Far Infrared Slab Lensing and Subwavelength Imaging in Crystal Quartz 
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We examine the possibility of using negative refraction stemming from the phonon response in 
an anisotropic crystal to create a simple slab lens with plane parallel sides, and show that imaging 
from such a lens should be possible at room temperature despite the effects of absorption that are 
inevitably present due to phonon damping. In particular, we consider the case of crystal quartz, 
a system for which experimental measurements consistent with all-angle negative refraction have 
already been demonstrated. Furthermore, we investigate the possibility of subwavelength imaging 
from such materials, and show that it should be possible for certain configurations. 

PACS numbers: 42.25.Bs, 42.25.Lc, 71.36.+C 



I. INTRODUCTION 



The idea of a slab lens stemming from negative re- 
fraction was described by Veselago as far back as 1968. 1 
Such a lens would have plane parallel sides, and an ob- 
ject O placed on one side of the lens would project a 
real image I\ within the slab followed by a second image 
I2 on the other side of it, in line with the well-known 
ray geometry shown in Fig. HJa). At the time, the con- 
cept, based on materials having both permittivity e and 
permeability fi simultaneously negative, was regarded as 
essentially hypothetical. Pendry's 2000 paper— brought 
it into the limelight, however, partly due to the realiza- 
tion that double negative materials (e < 0, fi < 0) were 
becoming a reality through metamaterial engineering, 3 
and partly as a result of the suggestion that the resulting 
lenses may have imaging possibilities beyond the tradi- 
tional diffraction limit. 

In considering refraction within a slab lens, it should 
be noted that it is the ray or Poynting vector direction 
(as opposed to the wavevector direction) that determines 
the focusing behavior. So long as negative refraction is 
defined in in this way, any slab of material that displays 
negative refraction for both positive and negative angles 
of incidence should display some degree of slab lensing 
behavior, regardless of the mechanism leading to nega- 
tive refraction. Thus, assuming the slab is sufficiently 
thick to create the intermediate image 1%, image forma- 
tion of the type shown in Fig. [IJa) (although not neces- 
sarily perfect imaging) should be possible within a cer- 
tain range of incident angles. One very simple way of 
achieving the necessary negative refraction is to make 
the slab from a nonmagnetic anisotropic medium two of 
whose principal axes have dielectric tensor components 
of opposing signs4^— In the correct configuration, such 
media induce negative refraction at all incident angles 9\ 
in the range —90° < 0\ < 90°, thus making them par- 
ticularly promising for slab lens construction^^ 3 - Suit- 
able materials may be constructed artificially^^ 3 - - — , 
but the required characteristics may also be present in 
certain natural anisotropic crystals around the phonon 



frequencies j ^ 10 ' 11 Recent experimental measurements 
on quartz have yielded results supporting the existence of 
all-angle negative refraction due to the phonon response 
at room temperature^ so it is natural to wish to consider 
how this type of material may perform as a slab lens. We 
consider such behavior in the present paper, as well as 
considering the possibility of using natural crystals such 
as quartz for subwavelength resolution. 



II. ALL-ANGLE NEGATIVE REFRACTION IN 
ANISOTROPIC DIELECTRIC MEDIA 

In order to understand slab lensing behavior in non- 
magnetic anisotropic media, we first consider how nega- 
tive refraction of a single ray may occur in the geometry 
shown Fig. [IJb). The slab is made of such a medium 
oriented with its principal axes along the cartesian axes 
x, y and z. xz is the plane of incidence and z is normal to 
the slab surface. If the angle of incidence is represented 
as 0i, the in-plane wavevector component k x is given by 



feo sin ( 



(1) 



where fco = w/c. In p-polarization (E field in the plane of 
incidence xz), the z component of the wavevector outside 
the slab is given by 



biz — k>3z — ^0 



(2) 



where the subscripts 1 and 3 represent regions to the left 
and to the right of the slab respectively. Inside the slab 
It2z is represented by 



(3) 



where e xx and e zz represent two of the principal compo- 
nents of the dielectric function of the anisotropic medium. 
The correct sign of k 2z is determined by the condition 
that power flow must be away from the interfaced 

Ray directions within each medium follow those of the 
Poynting vector S, so calculation of the vector S2 should 
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tion within the slab. A ray through the slab thus follows 
the Poynting vector directions shown in Fig. [TJb), with 
the wave behavior shown in Fig. QJc). Comparison of 
these two figures also confirms that the wavevectors re- 
main normal to the wavefronts in all layers. 



III. NEGATIVE REFRACTION DUE TO 
PHONONS IN NATURAL CRYSTALS 

One way of satisfying the condition e xx > 0, e zz < is 
to make use of the phonon response in natural anisotropic 
crystals! 7 ' 8 ! 10 " — In particular, we consider a uniaxial 
crystal whose dielectric tensor, expressed in relation to 
the crystal axes, takes the form 
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£ext 



(5) 



Here e cx t refers to the response along the extraordinary 
axis and e rd to the response along the ordinary axes. 
Around the phonon frequencies, these components may 
be written in the form 



(c) 




FIG. 1. (Color online) (a) Ray diagram of image formation 
from a slab lens, (b) Geometry considered in this paper show- 
ing wavevector and Poynting vector directions for a single in- 
cident ray. (c) Field profiles showing beam and wavefront 
directions for a single ray. 



yield the angle of refraction 6*2 within the slab. Its time- 
averaged value is given by (S2) = 1/2 Re (E x H*), lead- 
ing to a 62 value obtained from the expression 7 - 



tan 82 = 



2x1 
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Re(k x /e zz ) 
Re(k 2z /e xx )' 



(4) 



To a first approximation, we can ignore absorption, 
leaving both e xx and e zz real. In this case k,2 Z can be 
either real or imaginary, propagation into the slab occur- 
ring when it is real. We note in particular that this is 
the case when e xx > 0, e zz < 0, and it is straightforward 
to show that k 2z is positive under these conditions^ A 
simple comparison of Eqs. ([1]) and ^ thus shows that 
6\ and O2 have opposing signs, leading to negative refrac- 
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where oj is the frequency, e^ord and £oo,cxt are the 
high frequency dielectric constants, WTn,ord and WT„ i0X t 
are the frequencies of transverse optical (TO) phonons, 
^Ln.ord and WLn.ext are the frequencies of the longitudinal 
optical (LO) phonons, and 7T«,ord, 7Tn,ext, 7Ln,ord and 
7Ln,oxt are the appropriate damping parameters. In Fig. 
[5]we show the values of e 01 -d and e cx t for crystal quartz in 
the range 400 cm -1 to 600 cm -1 . The parameters used 
are based on those obtained by Gervais and Piriou J£ We 
have made adjustments to some of their values, however, 
in order to give a better fit to the experimental results 
presented later in this paper, as summarized in Table HI 

The dielectric function in the phonon region is, in gen- 
eral, complex, but it is reasonable, as a first approxima- 
tion, to simply to look at its real part in considering the 
refracting behavior. Thus, for negative refraction to take 
place, Re(e or d) and Re(£ oxt ) should have opposing signs. 
It is seen that, in the case of crystal quartz, Re(e rd) > 0, 
Re(e e xt) < in the frequency region between WL2,ord and 
WL2,ext (using the mode numbering of Table HJ whereas 
Re(eord) < 0, Re(e cxt ) > in the frequency region be- 
tween WT2,ord and WT2,cxt- The corresponding region of 
negative refraction depends on the crystal orientation. 

We first consider the extraordinary axis to be along z 
(i.e. normal to the crystal surface), so that e xx = e rd 
and e zz = e ox t. In this case we have the negatively 
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TABLE I. Comparison of phonon parameters of crystal quartz used in this work with those of Gervais and Piriou.— We have 
retained the high frequency dielectric constants used by these authors (eoo.ord = 2.356, Eoo.cxt = 2.383). 

Gervais and Piriou This work 









7Tn 




7Ln 




7Tn 




7Ln 


Symmetry 


n 


(cm *) 


/ —In 

(cm x ) 


(cm x ) 


(cm x ) 


(cm J ) 


(cm x ) 


(cm ) 


(cm ^ 




1 


393.5 


2.8 


402.0 


2.8 


393.5 


2.1 


403.0 


2.8 




2 


450.0 


4.5 


510.0 


4.1 


450.0 


4.5 


507.0 


3.5 


E 


3 


695.0 


13.0 


697.6 


13.0 


695.0^ 


13.0^ 


697.6^ 


13.01 


(ordinary) 


4 


797.0 


6.9 


810.0 


6.9 


797.6f 


6.9f 


sio.df 


6.9f 




5 


1065.0 


7.2 


1226.0 


12.5 


1065. qf 




1226. 


12.5f 




6 


1158.0 


9.3 


1155.0 


9.3 


1158. 0~ 


9.3~ 


1155.0* 


9.3~ 




1 


363.5 


4.8 


386.7 


4.8 


363.5^ 


4.8^ 


386.7 


7.0 


A 2 


2 


496.0 


5.2 


551.5 


5.8 


487.5~ 


4.0~ 


550.0 


3.2 


(extraordinary) 


3 


777.0 


6.7 


790.0 


6.7 


777.0^ 


6.7^ 


790.0^ 


6.7^ 




4 


1071.0 


6.8 


1229.0 


12.0 


1071 .(F 


6.8f 


1229. df 


12. of 






509.0 b 


14.0 b 


507.5 b 


14.0 b 











a The parameters associated with modes having frequencies outside the range of this study, or its immediate vicinity, have been left 

unchanged from those published by Gervais and Piriou. 
b This additional mode, used in modeling the experimental results of Gervais and Piriou, is generally considered to be an experimental 

artefact, 18 ! 19 and has not been included in our simulations. 
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measured using a Bruker Vertex 70 spectrometer, and 
theoretical simulations, obtained using standard trans- 
fer matrix techniques^ for various samples thicknesses I 
(see Fig. [T]) and incident angles Q\. 

We initially concentrate on the frequency marked as X 
in Fig. [31 at 531 cm -1 , since this frequency should yield 
negative refraction with relatively high transmission. We 
model the behavior of a finite beam passing through the 
slabs at this frequency by considering the incident beam 
to be gaussian, and represent it as a Fourier sum of plane 
waves: 
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FIG. 2. (Color online) (a) Real and (b) imaginary parts of 
the principal components of the dielectric function of quartz 
in the frequency range 400 cm -1 to 600 cm -1 . 



refracting condition Re(e xx ) > 0, Re(e zz ) < be- 
tween 507 cm -1 and 550 cm -1 (i.e. between 0JL2,ord and 
WL2,cxt), as shown in Fig. (3Ja). Negative refraction in this 
orientation has been studied in Ref. [H, which shows that 
significant transmission occurs in the corresponding fre- 
quency region. This is confirmed in transmission spectra 
shown in Fig. [31 which shows both experimental data, 



(8) 



In the case of a gaussian beam, ip(k x ) can be written^ 



ip(k x ) 



2 cos Oq^/tt 



exp 



g 2 (k x - k sin6> )^ 



4 cos 2 6 



(9) 



where 2g represents the beam width at its waist and 
6q represents the effective incident angle of the overall 
beam. In practice, we assume that all components of 
the gaussian beam are propagating in air (i.e. k\ z is 
real) j22 so we restrict the integral in Eq. © to the range 
—ko < k x < k Q . 

Using an incident beam of this form, it is possible to 
use standard multilayer optics techniques to calculate the 
magnetic field associated with each plane wave compo- 
nent at any point in the xz plane. Numerical integration 
then gives the overall H fields, and thus the associated 
E fields and Poynting vectorsi^ 2 ^ 

The resulting beam profiles for the various experimen- 
tal configurations represented in Fig. |3] are shown in 
Fig. 21 Here the incident beam, whose width is given 
by g — 100 /im, is assumed to be focused at the slab 
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FIG. 4. (Color online) Simulation of the intensity profile of 
a gaussian beam passing through a quartz slab in the con- 
figurations used in Fig. [3] at frequency X (531 cm" 1 ), (a) 
I = 25 (im, 0! = 30°; (b) I = 25 fim, (?i = 60°; (c) I = 50 yum, 
6>i = 30°; (d) I = 50 (im, f?i = 60°; (e) I = 75 fim, 0i = 30°; 
(f) I = 75 ^m, 6*i = 60°. 



FIG. 3. (Color online) (a) Simulated p-polarization angle of 
refraction in the case of quartz oriented with its extraordi- 
nary axis along z, for incident angles t?i of 30° and 60°; (b) 
transmission spectra through a crystal of thickness I = 25 fim 
at these angles of incidence; (c) transmission spectra through 
a crystal with Z = 50 fim; (d) transmission spectra through 
a crystal with I = 75 /»m. The solid and dashed lines in the 
transmission spectra represent experimental and simulated re- 
sults respectively. 



surface, at x = 0, z = 0. Negative refraction, seen as 
a displacement of the transmitted beam in the negative 
x direction in a manner similar to that shown in Fig. 
[IJc), occurs in each case. The displacement is naturally 
greater for thicker samples, but the transmission is lower, 
in line with the spectra shown in Fig. [31 In addition, the 



transmitted intensity is significantly reduced when the 
angle of incidence is increased. This is also observed in 
the experimental results. 



IV. SLAB LENSING IN NATURAL CRYSTALS 

We now consider how this kind of negative refraction 
may be used for slab lensing of the type shown in Fig. 
QJa). We take a source to be positioned at x = 0, z = 0, 
at a distance V to the left of the slab, i.e. the front 
surface of the slab is at z — V . As an approximation 
to a slit source, the amplitude of the incident H y field 
is assumed constant in the range —a/2 to a/2 at z — 
0, being zero at all other points in this plane. Thus a 
effectively represents a slit width. In a manner similar to 
that used to describe a gaussian beam, we represent the 
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FIG. 5. (Color online) Simulation of the intensity profile due 
to a slit source, of width a, at x — 0, z = 0. A slab lens 
of thickness I is positioned at a distance I' from the source. 
Calculations were made at frequency Y (537 cm -1 ), (a) I = 
25 fim, I' = 12.5 fim, a = 2.5 /jm; (b) { = 50 fim, I' = 25 /im, 
a — 5 /im; (c) i = 75 /jm, 2' = 37.5 /jm, a = 7.5 /jm. 



incident field to the right of the z = plane by Eq. ©, 
but ip(k x ) is now given by 

Mk.) = *B&g/*>. (10 ) 

The techniques used for calculating the overall fields and 
intensities in the xz plane are then the same as those 
used above for gaussian beam simulations. 

In the slab lens calculations, rather than the frequency 
X used in the simulation of transmission of gaussian 
beams, we use the slightly higher frequency marked as 
Y in Fig. |3[a), at 537 cm -1 . This is because, although 
the transmission is lower at this frequency, the angle of 
refraction O2 is (in magnitude) somewhat higher. The 
result of simulations at this frequency for the three slab 
thicknesses considered earlier are shown in Fig. [5] For 
each of these thicknesses we take V — 1/2 and a = l/W, 
so that, in a pure geometrical optics analysis that ignores 
diffraction limiting effects and absorption, the three fig- 
ures would be equivalent. Focusing of the two images I\ 
and I2 [see Fig. HJa)] clearly occurs in each case, but the 
image size does not simply scale with the overall dimen- 
sions of the system. 

There are several important considerations in the anal- 
ysis of the results shown in Fig. [5] Firstly, the slit 
width considered in the calculations are smaller than the 
diffraction limit in each case. In Fig. [SJa) its value is 



equivalent to a = 0.134A, in Fig. [5jb) to a — 0.267A, and 
in Fig. [5{c) to a — 0.403A. If the image size is diffrac- 
tion limited, it is therefore natural to expect the most 
pronounced increase in image size with respect to object 
size a in the situation shown in Fig. [5ja), where the 
object is smallest, and this is what is indeed observed. 
Secondly, the effects of absorption are not inconsiderable 
at the slab thicknesses considered, as observed in Figs. 
[3] and SJ This reduces the intensity of the image and, 
since the effect is larger for larger k x , may also adversely 
affect the image quality. Finally, even in a geometrical 
optics analysis, the image formation due to slab lenses 
with e xx > 0, e zz < is not perfect, and there are aber- 
rations associated with the higher incident angles. 7 

It is clear from the above that the images produced by 
slab lenses of the type considered here suffer from diffrac- 
tion limiting effects, and subwavelength details associ- 
ated with \k x \ > fcrjj corresponding to decaying evanes- 
cent fields in the air layers, are lost. Thus the lens is 
not functioning as a superlens of the type considered by 
Pendry^ In the Pendry lens, consisting of a slab with 
e = — 1, fj, = — 1, any decay in these fields is compensated 
by growing evanescent fields within the slab, so these de- 
tails are recovered. In the case of a nonmagnetic medium 
with e xx > 0, e zz < 0, in contrast, Eq. [3] shows that fez 
is real for all k x , i.e. there are no evanescent fields in 
the slab, either decaying or growing (although the prop- 
agating fields may suffer decay due to absorption). Thus, 
this type of lens can never operate in the same way as a 
Pendry lens in achieving subwavelength resolution. 



V. SUBWAVELENGTH IMAGING 
POSSIBILITIES 

As discussed above, restoration of evanescent waves is 
not possible for this type of lens, since evanescent waves 
are not present within the slab. However, the absence of 
evanescent waves may be used to advantage if the object 
is placed at the left hand surface of the slab, i.e. V = 
0. We can then consider an image, formed solely from 
propagating waves, at the right hand surface of the slab. 

The Pendry lens relies on the principle that the overall 
optical path from object to image is zero for all k x . In 
fact, to fully restore the object details, it is sufficient 
that the optical path should be independent of k x , i.e. 
the overall phase change should be constant. In the case 
of an object at the left hand side of the lens and an 
image at the right hand side of it, this (ignoring phase 
changes across the the interfaces) simply amounts to ki z 
being independent of k x . As seen from Eq. [31 this occurs 
when e xx > 0, l/e zz — > 0, so under these conditions 
subwavelength imaging should occur. 

It is noticeable that the condition l/e zz — > merely 
requires that the amplitude of e zz be large, without any 
restriction on its sign. In fact, it is not even required 
to be real, so a large imaginary e zz satisfies the condi- 
tion. From Fig. [21 we see that, within the range investi- 
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FIG. 6. (Color online) Oblique incidence p-polarized trans- 
mission spectra through various thicknesses of quartz crystals 
having their extraordinary axes directed along x. The angle of 
incidence is 9\ is 30 in each case. The solid and dashed lines 
represent experimental and simulated results respectively. 



gated, the combined condition e xx > 0, l/e zz — !• does 
not occur in quartz if the extraordinary axis lies along z 
(£xx = £ord, £zz = £oxt), but occurs at the TO phonon fre- 
quency wx2,ord (450 cm -1 ) if the extraordinary axis lies 
along x (e xx = e cxt , s zz = e Q rd), since lm(e zz ) becomes 
large at this frequency. This is therefore the geometry 
we use in our discussion of subwavelength imaging. 

Since we are considering transmission across the slab 
in p-polarization at a resonance frequency, we should 
first check that there is no absorption associated with 
this resonance. In Fig. [5] we show both experimental 
and theoretical p-polarized transmission spectra in the 
required geometry at oblique incidence. It is clearly seen 
that there is no absorption dip at WT2,ord; even though 
the z-component of the incident E field is nonzero in p- 
polarization. We can interpret this in the following way. 
Boundary conditions dictate that D z should be contin- 
uous across the interface, so a large \e zz \ implies that 
E z — >• in the slab, and the TO mode is not excited. 
There is some absorption in this region, as can be ob- 
served from the decreasing transmission with increasing 
slab thickness, but this is essentially due to the proximity 
of the s-polarized TO phonon at WT2,ext- 

Figure [7] shows the amplitude and phase of the H y field 
transmitted through slabs of quartz of the three studied 
thicknesses as a function of normalized in-plane wavevec- 
tor k x /ko. In the ideal case, both amplitude and phase 
would be constant for all k x . In practice, there are no- 
ticeable deviations from this ideal behavior. 

Figure Ul. a ) shows that there is a considerable drop off 
in amplitude at large \k x \. We can regard this as due 
to two separate effects, transmission efficiency across the 
two interfaces at either side off the slab and absorption 
within the slab. The first effect gives the basic shape of 
the curves and the second effect accounts for the separa- 
tion of the three curves representing the three different 
thicknesses. It is clear that the first effect dominates 
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FIG. 7. (Color online) (a) Amplitude and (b) phase of the 
H v field transmitted through various thicknesses of quartz 
crystals having their extraordinary axes directed along x, as 
a function of k x /ko, in p-polarization. 



the curves. Figure [7fb) shows that there is some phase 
change with k x , but the overall variation for a partic- 
ular slab thickness is around ir/2 in the range shown. 
Similarly to the amplitude curves discussed above, the 
basic shape is associated with phase changes on transmis- 
sion across the interfaces. Phase changes associated with 
transmission within the slab simply give a vertical shift to 
this basic shape, since these phase changes are almost in- 
dependent of k x . From the above, we see that the main 
restrictions to the required subwavelength imaging be- 
havior are likely to be associated with transmission across 
the interfaces. A number of studies of the use of metallic 
layered structures to achieve the required anisotropic di- 
electric tensor, have also discussed this phenomenon^r— 
Of particular importance is the result suggesting that 
use of a slab thickness equal to an exact number of half- 
wavelengths (i.e. k2 Z l = run where m is an integer), 
thus assuring constructive interference from Fabry-Perot 
fringes, should overcome these restrictions ! 28 ' 29 A special 
case of this, equivalent to choosing m = 0, is possible if 
£xx = 0. In the present work, we can see from Fig. [5] 
that some weak interference fringes are observed in the 
transmission spectra in the case of the 25 ^m-thick sam- 
ple, but that they are essentially absent in the case of the 
thicker samples. Thus we believe that, for the range of 
thicknesses considered in this work, the Fabry-Perot con- 
dition is not of crucial importance since the higher order 
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partial rays are absorbed by the slab. In fact, the 25 jam- 
thick sample is close to satisfying the Fabry-Perot condi- 
tion with m = 6 (an exact calculation gives m — 6.13), 
but a small change in the slab thickness does not appear 
to have much effect on the results. 

We now turn to simulations of subwavelength imaging 
itself. We consider a two-slit source in which the mag- 
netic field of the incident beam is assumed to be constant 
across the width of each slit, as before. For slits of width 
a separated by a distance d, this amounts to setting ip(k x ) 
to 
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2 sin(fca;a/2) cos(k x d/2) 
7rfc x 



(11) 



The slits are placed at the front surface of the slab 
(V = 0). Figure E] shows the resulting intensity distri- 
butions in the case of the 25 /im-thick slab. Figure Eta) 
shows results for a slit separation of d = 7 /jm (0.32 A) 
and slit widths a — 2.5 /im (0.11 A). The intensities from 
the two slits are well resolved within the slab, with sig- 
nificant loss of intensity with propagation through the 
slab. The solid (green) curve in Fig. [9ja) shows the in- 
tensity distribution passed through the slab. The images 
from the two slits are still easily resolved. When the slit 
separation d is reduced to d = 5 fim (0.23 A), we have 
found that the images are better resolved if we also re- 
duce the slit widths. We therefore show the intensity 
distribution for a slit separation of d = 5 /im and slit 
widths a = 1.5 fj,m (0.07 A) in Fig. Mb). The intensities 
from the two slits are still well resolved within the slab, 
but are considerably reduced due to the narrowing of the 
slits. The intensity distribution passed through the slab 
is shown as the dashed (red) curve in Fig. H])(a), and the 
two peaks are once more resolved. 

Figures^ (b) and[H])(c) show the intensity distributions 
passed through slabs of thicknesses I = 50 /im and I = 
75 /im respectively, using the same slit width/separation 
combinations as for the thinner slab. It is seen that when 
I = 50 /urn the images of the slits are still resolvable for 
d = 7 /im, but not for d = 5 /im. When I = 75 /urn, some 
structure remains in the d = 7 /im case, but not in the 
case of d = 5 /im. 

We thus see that subwavelength imaging should occur 
even for relatively thick slabs of quartz, corresponding 
to a few free-space wavelengths, albeit with considerable 
loss of intensity. 



VI. CONCLUSIONS 

The above results confirm that simple anisotropic crys- 
tals, such as quartz, should function as slab lenses as 
well as achieving images with subwavelength resolution. 
We have restricted our simulations to the slab thick- 
nesses used in our experimental spectral investigations, 
and such thicknesses are easily obtainable (the quartz 
crystals used in the spectroscopic measurements were ob- 
tained commercially from Boston Piezo-Optics). 




FIG. 8. (Color online) Simulation of the intensity profile due 
to a two-slit source, of width a and separation d, at the sur- 
face of a 25 /im-thick slab of quartz with its extraordinary 
axis along x. Calculations were made at frequency u>T2,ord 
(450 cm" 1 ), (a) a = 2.5 /im, d = 7 /im, (a) a — 1.5 /im, 
d — 5 /im. 



The slabs used for our subwavelength studies are fairly 
thick in relation to those in most studies based on mul- 
tilayer structures Better resolution should be pos- 
sible with thinner slabs, and it should also be possible 
to take advantage of Fabry-Perot interference in such 
casesj 28 ' 29 Nevertheless, in practice, the behavior at the 
interfaces may be drastically affected by the source and 
detector configuration if they are close to the surfaces, 
and a plane- wave analysis, although correct within the 
slab, may not give an accurate indication of the interface 
behavior. 

In this study we have only considered crystal quartz as 
the slab medium, but there are a number of anisotropic 
crystals that may be suitable. Several factors may be 
important in choosing suitable materials. Obviously 
different materials will be appropriate for different fre- 
quency ranges, and phonon resonances must be suffici- 
iently strong and well separated. Absorption clearly 
plays a vital role in the image formation, so low damp- 
ing is important. For subwavelength imaging using the 
£xx > 0, l/£zz — > criterion at z-polarized TO fre- 
quences, it would be useful to have negligible lm(e xx ), 
so there should ideally be no any a;-polarized phonons 
close to the frequency of interest. Dvorak and KuzelS. 
discuss the case of Hg2l2 in the context of negative re- 
fraction (rather than imaging behavior). The damping 
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FIG. 9. (Color online) Intensity profile transmitted across a 
slab of quartz with its extraordinary axis along x, due to a 
two-slit source of slit width a and separation d. Calculations 
were made at frequency UT2,ord (450 cm" 1 ). The intensity 
scale is normalized with respect to a plane wave whose mag- 
netic field amplitude is that of the incident field in the slits, 
(a) Slab thickness I — 25 ^m, (b) I = 50 fim, (c) I — 75 /im. 



parameters for this material are somewhat larger than 
those for quartz, but its optical phonon frequencies along 
the principal axes are well separated from one another, 
so may be suited to subwavelength imaging applications. 
Dumelow et al£ consider slab lenses from triglycine sul- 
fate (TGS), which, at low temperature, has both very 
low damping and well-separated phonon frequencies. In 
addition to the slab lensing properties discussed in the 
paper, it is likely to give very good subwavelength imag- 
ing. The disadvantage is, of course, the necessity of low 
temperature. Another material considered for this type 
of study is MgF2, which has properties somewhat similar 
to quartz iiS 

In this type of study and its in subsequent applica- 
tions it would be useful to observe single-frequency be- 
havior experimentally. There have been recent reports 
of quantum-cascade lasers operating in the frequencies 
discussed in this paper; 30 ! 31 which at present are at 
the lower frequency limit of this technology in the in- 
frared region. In the immediate future, spectroscopic 
measurements, with the aid of stops, slits or gratings, 
may be an easier option for investigating imaging effects 
in quartz. However, optical phonons, in general, span 
a wide frequency range, and investigations of imaging 
properties using monochromatic radiation may be more 
straightforward in other materials. The higher frequency 
modes of calcite^ for example, should be far more 
easily accessible than those of quartz using quantum- 
cascade lasers. At considerably lower frequencies, below 
~ 150 cm -1 , quantum-cascade lasers may again be con- 
sidered as possible sources^ along with other devices 
such as backward- wave oscillators. 34 Imaging using crys- 
tals such as Hg2l2 an TGS, whose phonons fall in this 
frequency range, should therefore be possible with the 
aid of such sources. Overall there appear to be ample 
possibilities for using anisotropic crystals in this way. 
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